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O !• Introduction 

(N ■ 

The study of general Hamiltonian structures for systems of hydrodynamic type, i.e., 
^ . systems of the form 

JV 

«* = J>i(u)«;, i = i,...,N, (i) 

k=l 

was initiated by Dubrovin and Novikov [2] and continued by Mokhov and Ferapontov [5] . In 
the present paper, we prove a theorem on the existence of three Hamiltonian structures for a 
diagonalizable Hamiltonian hydrodynamic type system (1) having two physical symmetries, 
with respect to the Galilean transformations and to scalings, and possesses some additional 
properties, namely, the metric of the Hamiltonian structure has the Egorov property and 
the matrix ((r l — r k )(3 ik ) (see Sec. 4 below) is semisimple. There are rather many known 
physically meaningful systems of hydrodynamic type which have this form, including av- 
eraged equations on iV-phase solutions of the Korteweg-de Vries equation (the Whitham 
equations) and the nonlinear Schrodinger equation as well as the dispersionless limits of 
the vector nonlinear Schrodinger equation and the vector long-short resonance equation. For 
the Whitham equations, local multi-Hamiltonian structures were described in [15, 1] by the 
algebro-geometric method. In the present paper, we construct third Hamiltonian structures 
for all above-mentioned systems by a common differential-geometric method. It turns out 
that to find out whether there exist second and third Hamiltonian structures and determine 
their type (local, nonlocal with a constant curvature metric, or general nonlocal), it suffices 
to know the homogeneity degrees of the annihilators of the first Hamiltonian structure. These 
^ ■ degrees are usually a priori known. We also prove a simple criterion for the Egorov property 
of a metric (Theorem 1). Theorems on tri-Hamiltonian structures are stated and proved for 
the general case in Sees. 2-4 and are then applied to the above-mentioned examples in the 
subsequent sections. 

The theory considered here can be applied not only to integrable diagonalizable systems 
of hydrodynamic type but also to the theory of Frobenius manifolds and conformal topo- 
logical field theory for the solution of the Witten-Dijkgraaf-Verlinde-Verlinde associativity 
equations [18]. 

Let us briefly state some results needed in the sequel (see [13]). 

The original diagonalizable hydrodynamic type system (1) has Riemann invariants r l = 
r l (u), i = 1, . . . , N, i.e., variables in which the velocity matrix (v l k ) of system (1) is diagonal 
(unless otherwise stated, summation over repeated indices is not assumed): 

ri = v\r)rl z = l,...,iV. (2) 
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In the following, we assume that v l ^ v k for i ^ k and that system (2) does not split in the 
sense that diV k 7^ 0, i 7^ k. Proceeding from v l (r), we find the Lame coefficients Hi{u) as 
some solutions of the overdetermined system 

d k \nH t = Tl k = -^-, i^k. (3) 

Here d k = d/dr k and the T\ k are the Christoffel symbols of the Levi-Civita connection 
corresponding to the metric entering the Hamiltonian operator of the system (see (8) 
and (9) below). We also find the rotation coefficients (see [13, 17]) 

0* = i*k. (4) 

The compatibility conditions for systems (3) or (4) have the form of the semi- Hamiltonian 
property [12, 13] 

diP jk = Prf* ^ d 3 = d k , i^j^k. (5) 

Moreover, the diagonal coefficients of system (2) can be represented in the form v l (r) = 
Hi/ Hi, where the Hi are the Lame coefficients and the Hi are some solutions of the system 

d t H k = i5 ik H { , 1 ± k. (6) 

A local Hamiltonian structure of Dubrovin-Novikov type exists for a semi-Hamiltonian sys- 
tem (2) if and only if 

Wit + 9kP ki + E Pr*Pn* = 0, i^k. (7) 

m^i, k 

In this case, the hydrodynamic type system (1) or (2) can be rewritten in the Hamiltonian 
form 

A 6H ' 



u\ = K, H} = ^A^— = (8) 
3=1 k=i 

with Hamiltonian operator 

^ = g^u)^--Y,9 ts {n)T{ k {u)u k x , (9) 

s,k 

where g l i = is a nondegenerate zero curvature metric and the connection T J sk is symmetric 
and compatible with the metric. (Here H = J h(u) dx is a hydrodynamic type Hamiltonian.) 
In the Riemann invariants, the metric is diagonal, i.e., gu = Hf (g n gu = 1), and the vanishing 
of the Riemann tensor is equivalent to (7). Thus, with each local Hamiltonian structure of 
system (2), a system of orthogonal curvilinear coordinates r % in a flat iV-dimensional pseudo- 
Euclidean space is associated. Note that some of the functions Hi and f3 ik are pure imaginary 
in the case of pseudo-Riemannian diagonal metrics. (It is metrics of this type that occur in 
applications, as is elucidated by Theorem 7 below.) 

Each coefficient Hi is determined by (3) to within a multiplication by a function of one 
variable, Hi = ^^(r^Hf^. In this case, condition (7) acquires the form 

\ ^f k ] + \ ^£ X) + M/4 fix) + »MF + E v^JS^S = o- (io) 

m^i, k 



Since system (10) is linear in fj, iy we can readily see that any two Hamiltonian structures (9) 
associated with a given diagonal hydrodynamic type system with nondegenerate metric gij 
are automatically compatible. As was shown in [20], system (10) for /ij(r l ) can have at most 
N+ 1 linearly independent solutions. In the general case, it has no solutions at all for a given 
semi-Hamiltonian system. 

The hydrodynamic type system (1) can also have the nonlocal Hamiltonian structure (8) 
with a nonlocal Hamiltonian operator [11] 

A« = g ij d/dx - £ 9 ls nA + £ s^w^u^d/dxy'w^u™. (11) 

s,k a,/3,n,m 

Here the wlf a \u) are the matrix coefficients of the hydrodynamic type flows 

N 

M U = E w i (a) ( M K (12) 

commuting with (1). In the Riemann invariants, the coefficients w^ a ' are diagonal, namely, 
w h{°<) _ Hi, where h\°^ is a set of solutions of system (6) such that 

d^k + d k (3 kl + £ PmiPmk = E ^^H^, i ± k, (13) 

£ a /3 = £/3a = const. Thus, the metric g^iu) in the above nonlocal Hamiltonian operator is 
not flat. If Eq. (13) contains a single term HiH k formed of the Lame coefficients if; = v /(fe 
of the given metric, we obtain a constant curvature metric [5]. 

2. The Egorov Systems of Hydrodynamic Type 

As is shown in [12,13], a diagonal semi-Hamiltonian system (2) has infinitely many 
commuting hydrodynamic type flows 

rl = w\ry x , i = l,...,N, (14) 

whose coefficients w % (r) can be found as solutions of the consistent overdetermined system 
of first-order linear partial differential equations 

d k w l = T\ k {w k -w l ), i^k. (15) 

System (2) also has infinitely many hydrodynamic conservation laws dp(r)/dt = dq{r)/dx, 
where the conservation law densities p(r) can be found from the consistent system 

d l d kP = T\ k d l p + T k ik d k p, i^k. (16) 

The substitution w l = Hi/Hi reduces system (15) to the form (6), and the substitution 
dip = i>iHi reduces (16) to the system 

dki>i = Pik^k, i^ k , (17) 

which is the adjoint of (6). 

Definition 1. An orthogonal coordinate system associated with the diagonal metric 
ga = Hf is called a The Egorov coordinate system if 

Pik = Pki- (18) 

In this case, the metric is potential, i.e., ga = dia(r) for some function a(r). 
The corresponding hydrodynamic type systems will also be called The Egorov systems. 
Condition (18) is not invariant with respect to the natural transformation r l — > (p l (r l ) of 
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Riemann invariants or the transformation Hi — > 1 ^ 2 (r l )H i . An invariant condition for a 
nonsplitting system (f3 ik ^ 0) with N > 3 (see [3, 17]) is given by the relation 

PikPkjPji = PijPjkPki, i^j^k. (19) 

For a Egorov system of hydro dynamic type, the relations d(p = HiHi hold and systems (6) 
and (17) coincide. Integrable systems of hydrodynamic type have infinitely many indepen- 
dent conservation laws d t pk(u) = d x qk{u), k = 1,2, ... , and hence can be rewritten in the 
conservative form 

d t a a = d x q a (a), a = l,...,N, (20) 

where the a a = a a {u) are conservation law densities and the q a (a) are the corresponding 
flows. 

Let us prove the following criterion for the Egorov property of system [8]: 

Theorem 1. A diagonalizable nonsplitting semi-Hamiltonian system (1) has a pair of 
conservation laws of the form 

a t = b x , b t = c x (21) 
if and only if it is a Egorov system. In this case, in the Riemann invariants, the relations 
dia = Hf , djb = HiHi, and diC = H i hold, where Hi = v l (r)Hi. 

Proof. Condition (21) can be rewritten as 

d t a = ^Hi, d t b = ^Hi, (22) 
dib^Hi, «9,c = CT, (23) 

where the ipi an d V'i are solutions of system (17), whence it follows that ipiHi = ipiHi. 
Differentiating tp i = ip^i/Hi with respect to the Riemann invariant r k (i ^ k) and dividing 
the result by V{ — v k = Hi/ Hi — H^/ H k , we obtain 

^/?* = ^Ak (24) 
ilk tli 

from (17), which implies (19) provided that the number N of equations is greater than 2. 
Note that, by virtue of (3) and (4), the nonsplitting condition [diV k ^ 0, i ^ k) implies 
that P ki 7^ for % ^ k and by virtue of (24), it guarantees that t/fj ^ and dia ^ 
for any i. To prove the Egorov property for the case in which system (1) contains only 
two equations, we set qi = ipi/Hi = ifj i /H i . It follows from (23) that di(qkHf) = d k {qiHf) 
and di{qkH k Hk) = dk{qiHiH,j) , whence, using (24), we arrive at the relations Hfd k qi = 
H\diq k and HiHid^qi = HkH^diq^. Dividing one of these expressions by the other (under 
the assumption that d^qi ^ 0, i ^ k), we obtain Vi = Vk, which contradicts the assumption 
that Vi 7^ Vk for i ^ k. Consequently, d k qi = 0. Since the Riemann invariants are defined only 
to within the transformation r l — > (p l (r l ), we can always make the functions qi{r l ) = ipi/Hi 
in (24) identically equal to unity. Thus, the metric is potential and we have dia = Hf, 
dfi = HiHi, and <9jC = H { . 

Conversely, if the metric is potential, i.e., Hf = dia, then conditions (18) hold. In this 
case, we have a t = dia ■ r\ = Hfr\ = Hfv l (r)r l x = HiH(r x = b x and b t = dfi ■ r\ = HiH^l = 

H i H i v i {r)r i x = I? i ri = c x . | 

Note that the condition that the metric is flat has not been used in the proof of the theo- 
rem. Egorov himself studied potential metrics with the additional zero curvature condition. 
In this paper, potential metrics of arbitrary curvature will be referred to as Egorov metrics. 
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Since all commuting flows have the same metric (see (3) and (15)), we arrive at the following 
result. 

Corollary 1. An arbitrary commuting flow (14) of a Egorov semi-Hamiltonian system 
of hydrodynamic type has a pair of conservation laws of the form (21) d y a = d x h, d y h = d x g. 

It is obvious that if the hydrodynamic type system (2) has the local Hamiltonian struc- 
ture (8), then the relationship w l k (u) = V l VkP between the solutions of Eqs. (16) for conser- 
vation law densities and the solutions of Eqs. (15) for commuting flows holds. Writing the op- 
eration of covariant differentiation V* in full and substituting the expressions w l (r) = Hi/Hi 
and dip = ip^Hi (in Riemann invariants), we obtain a relation between the solutions of prob- 
lem (6) and the adjoint problem (17), 

H^d^ + ^Pmi^m- (25) 

If a hydrodynamic type system has a second local Hamiltonian structure, then, as was 
already noted, the diagonal coefficients g l fo of its metric can differ only in factors depending 
on the corresponding Riemann invariants, g%^ = ^i{r l )g l Ay Let us make a transformation 
r % — > (p l (r l ) of the Riemann invariants that gives l^i(r l ) — > r l (under the assumption that 
/i-(r l ) 7^ 0), i.e., gjL. = r % g l Ay By analogy with (25), for the connection generated by the 
second flat metric g l fa we obtain 

Hi = \il>i + + J2 r m (3 mi *p m . (26) 

If g % fa is an Egorov metric, then formulas (25) and (26) imply that 

Hi = 6^, ^ = (R+ 1/2)^, (27) 

where 5 = ^ is a translation operator and R — r k dk is a scaling operator (see [10]). 

Let us summarize the preceding. For an appropriate choice of Riemann invariants, if 
an Egorov hydrodynamic type system has one local Hamiltonian structure, then 5j3 ik = 0; 
in other words, the rotation coefficients of orthogonal coordinate systems depend only on 
the differences of Riemann invariants, f3 ik = [3 ik (r m — r n ). If the Egorov hydrodynamic 
type system has also a second local Hamiltonian structure with metric g l fo = r l g™ x y then 

Eqs. (10) and (5) imply the equivalent homogeneity condition R(3 ik = —f3 ik ; thus, the rota- 
tion coefficients of orthogonal coordinate systems are homogeneous functions of order — 1. 
As is shown in [13], the Egorov property of the metric for a diagonalizable Hamiltonian 
system essentially expresses the fact that the system is invariant with respect to the group 
of Galilean transformations (x, t) — > (x — V - t,t). The homogeneity of rotation coefficients in 
some well-known examples corresponds to the natural physical condition of invariance with 
respect to the choice of units of measurement, which results in the homogeneity of Vi(r). 
Thus, the presence of two local Hamiltonian structures for physical examples of hydrody- 
namic type systems is a manifestation of the simplest fundamental invariance conditions. We 
have already noted that the corresponding Hamiltonian structures prove to be compatible 
automatically. 

In the following, for brevity we say "homogeneous rotation coefficients" , implying that the 
degree of homogeneity is always —1, Rf3 ik = —f3 ik . The corresponding curvilinear orthogonal 
coordinate systems and diagonal hydrodynamic type systems (2) will also be simply said to 
be homogeneous. 
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3. Annihilators of the First Local Hamiltonian Structure and the Signature of 
its Metric for the Homogeneous Egorov Systems 

Annihilators of a local Hamiltonian structure are denned as conservation law densities 
a a (u) determined by the vanishing condition for commuting flows: 

= V l V fc a a , a = l,...,N. (28) 

The overdetermined system (28) has an iV-dimensional solution space, which specifies flat 
coordinates for the metric gu = Hf of the Hamiltonian structure. Introducing ip^ = did 01 /Hi 
in the Riemann invariants, we rewrite Eq. (28) in the form 

dk^ ] ={3 lk ^\ = % (a) + ^/5 m ^ ) , a = l,...,N. (29) 

The following two theorems were proved in [14] and later presented in [18]. 

Theorem 2. The homogeneous solutions ijjf*^ of system (29) for an Egorov orthogonal 
coordinate system whose rotation coefficients are homogeneous and depend only on the dif- 
ference of the Riemann invariants are eigenvectors of the skew-symmetric matrix B with 
entries = {r l — r J )/3^- ; which has constant eigenvalues. 

Proof. We have 

The constancy of the eigenvalues c a (r) of the matrix B can be proved as follows: di(Rip^) = 

R&fP) + MTJ = R(P ki 4 a) ) + Mi a) = R(PkiW ] + 0*%i a) ) + MP = -Mi a) + 

Ca^Pki + Pki&i^ = C a^\ a) (3kn i- e -> ^(ca-0fe Q) ) = di(c a )ijj k a) + c a (3 w ^ a) = c a ^°° '(3 ki , whence 
di(c a ) = 0. | 

The expressions ip^ will also be referred to as annihilators. 

Since the homogeneity degrees of any conservation law densities, including annihilators, 
are real numbers in the case of homogeneous real hyperbolic systems (2) considered here, 
it follows that the eigenvalues of the skew-symmetric matrix B are also real. This can be 
accounted for by the fact that for this class of hydrodynamic type systems the metric is 
necessarily pseudo- Euclidean (see Theorem 7 below), i.e., some Lame coefficients Hi = ^fg~^ 
and rotation coefficients /3 ik are pure imaginary. Thus, the matrix B has both real and pure 
imaginary entries, and one can no longer claim that there is a complete set of eigenvectors. 
For example, if N — 3, then for an Egorov homogeneous metric of signature 1 one can 
choose coefficients (3^ such that (r 1 — r 2 ) 2 l3\ 2 + (r 1 — r 3 ) 2 (3l 3 + (r 2 — r 3 ) 2 (3l 3 = 0. (It is easily 
seen that this condition is compatible with Eqs. (30) for f3 ik) given below.) In this situation, 
the matrix B has one eigenvector and two root vectors with eigenvalue 0. Note that the 
set of The Egorov orthogonal coordinate systems of given dimension with homogeneous 
rotation coefficients depending only on the differences of Riemann invariants is a finite- 
parameter family, and the description of such systems can be reduced to that of solutions of 
the consistent Pfaff system 



9iP jk = PjiPik, 5P ik = 0, Rf3 ik = -f3 ik) i^j^k, 
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(30) 



whose solution can be specified by setting an arbitrary symmetric matrix ((3 ik ) (where the 
f3 u are undefined) at the initial point (V( ))- 

However, all physical examples of diagonalizable homogeneous Egorov systems known 
to the authors, as well as The Egorov orthogonal coordinate systems arising in conformal 
topological field theory, for which the rotation coefficients are homogeneous and depend 
only on the differences of Riemann invariants, result in matrices B with complete sets of 
eigenvectors. Moreover, as a rule, the homogeneity degrees of annihilators of the first local 
Poisson bracket are a priori known for physical systems. Hence we assume throughout the 
following that the semisimplicity condition holds, namely, the matrix B has a complete set 
of eigenvectors at some point (r^). As was shown in [18], the description of semisimple 
solutions of system (30) can be reduced to the solution of the sixth Painleve equation. 

Theorem 3. For a semisimple matrix B , all annihilators corresponding to a flat 
homogeneous Egorov coordinate system can be taken in the form of homogeneous functions, 

->(a) ^->(a) -'(c) 

Bi > =Ri > =c a i - , a = l,...,N, (31) 

where the c a (the homogeneity degrees of the annihilators) are simultaneously the eigenvalues 
of the matrix B . 

Proof. For a symmetric matrix (/3 ik ), i.e., for The Egorov metrics, the second relation 
in (29) becomes 5^ = d^ a) + Yn^Mft = d^f + £ m ^ d m ^ a) = 0. Thus, the 
(possibly nonhomogeneous) functions ijjf^ satisfy the system 

S4 a) = 0, d^ k a) = P lk 4 a \ i*k. (32) 

One can readily verify that this system is consistent and the initial data for it are given by 
iV constants, namely, the values of the function $ at some point (rfo). Reproducing the 

computations made in the proof of the preceding theorem, we see that di(Rip k — ct/j k ) = 
l3 ik (R>p i — cipi) and S(Rip k — cip k ) = for each solution ip i and each c = const. Thus, 
ip i = Ripi — cipi also satisfies system (32). Consequently, equipping system (32) with the 
initial data ip ( r (o)) forming a basis of eigenvectors of the matrix B at this point (rLj, we 

see that the ip i = Rip i — c a tlj { are also solutions of system (32) that vanish at the point 

~(«) 

(r^). Thus, ip { = 0, as desired. | 

Definition 2. The momentum density P for a system of hydrodynamic type with local 
Hamiltonian structure is understood as the quadratic expression 

P = \Y,9^a a a\ (33) 

where the nonsingular matrix g a p is the (constant) metric of system (2) in flat annihilator 
coordinates a a . 

It is obvious (see [2]) that in the annihilator coordinates a a the hydrodynamic type 
system (1) has the form 

v (3=1 7 
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where h(u) is the Hamiltonian density. The momentum density gives an integral of an arbi- 
trary hydrodynamic type Hamiltonian system (34) with given metric g alS and generates the 
trivial commuting flow u\ — u l x . 

Note the following assertion (proved independently in several papers by different authors). 

Theorem 4. // a {not necessarily diagonalizable) hydrodynamic type system (1) has 
N + 1 linearly independent hydrodynamic conservation laws and one of their densities can be 
expressed quadratically via the others {see (33)), then the system has the local Hamiltonian 
structure (34) with constant metric g a p. 

Proof. Let d t a a = d x q a and d t P = d t [\ J2 a p g a pa a a^\ = d x Q(a). It follows that 
dQ = g a( 3a a dqP{a) = d[g a pa a q p ] - ^ g a pq P da a . 

a, 8 a,3 a,B 

Consequently, Ylsdasq 13 = dh/da a for some function h(u), and hence system (20) acquires 
the form (34). | 

The following theorem, generalizing Theorem 4 to the case of a constant curvature metric, 
was stated in [7] and proved in [25]. 

Theorem 5. Suppose that a hydrodynamic type system represented in the form of con- 
servation laws d t c a = d x b a has an additional conservation law density p quadratically related 
with the densities c a : 



(35) 



a,/3 



where g a p is a constant nonsingular matrix. Then in the field variables c a the system has 
the nonlocal Hamiltonian structure 



d x 



N 



dh 



Y.^-ec^)—p+ec-h 



L /3=l 



(36) 



associated with the metric g a/3 = g a/3 — ec a d 3 of constant curvature e, where (g al3 ) is the 
inverse matrix of {g a (i) ■ 

Proof. It follows from the formulas d t c a = d x b a and p t = d x Q(c) and from (35) that 



E c /3 
9a(S- db a . 
l — ep 

a,/3 



Pt = d x 



E 



ep 

and cf = d x 



By setting q a = c a /(l — ep), we find that J^ggasb^ = ds/dq a with some potential s. Hence 
we see that 



ds 
dq~? 



Since 



ds 



= {l-ep) 



ds 



--E 



0,1 



7 ds - t 
^c r - — g a RC 



we readily obtain (36) with h — (1 — ep)s. | 

Needless to say, the shift p — > p + const makes relation (35) purely quadratic. However, 
the form given in (35) clarifies the passage to the limit as e — > 0. 



Theorem 6. The following relations hold for the components of an arbitrary (not nec- 
essarily Egorov or homogeneous) flat metric in flat annihilator coordinates: 



N 



a,P i=l 

For the special case of an Egorov metric, this result was obtained in [19]. 
Proof. The momentum density P (see (33)) satisfies the relation 

diP = * i H i = Y t 9o0a P 4 a) H i => * i = Y l 9apa f> 4 a) - 

a,P ot,(3 

Differentiating it with respect to r k and using (17), we obtain - gap^i*^^ 
Taking account of the fact that P generates the flow u\ = u l x and differentiating with 
respect to r l , we obtain Yl a g gaP^t*^^ = 1> which precisely gives the first formula in (37). 



(37) 



0. 



Multiplying it by ip^ and performing summation over k, we arrive at the relation 



E ( E s^-Vf ) *f = E ^ ( E « E *iW 

fc=l ^ a,0 ' a=l ^ /3=1 k=l 

7(a) 

Since the vectors ip are linearly independent, it follows that 

N / N 



1> 



(7) 



EmE^wW, 

a=i ^ fc=i ' 



P 

i.e., we have arrived at the second formula in (37). | 

Let us return to the case of an Egorov orthogonal coordinate system whose rotation 
coefficients are homogeneous and depend only on the differences of Riemann invariants. 
Then 

N N N N 

Rg<* = o = Rj24 a W = wf } ) + wS Q) ) = ^ + <*) I>S Q Vf. 

i=l i=l 1=1 1=1 

It follows that either c a + Cp ^ 0, and then the corresponding component of the matrix 
(g a/3 ) is zero, or c a + cp = 0, and then the corresponding component of the matrix (g alS ) 
is nonzero. (Some of the components can also be zero, but there must be sufficiently many 
nonzero components to ensure that the metric is nondegenerate.) Hence if the semisimple 
matrix B has at most one eigenvector with zero eigenvalue, then, after an appropriate 
renumbering of the flat coordinates corresponding to the homogeneous functions the 
flat metric acquires the block antidiagonal form 

/ A \ 



a/3 







* * 







-.otP 







* 





±1 






* 







(38) 



N = 2n 



N = 2n + 1 



Thus we have proved the following assertion. 
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Theorem 7. // the multiplicity of the zero eigenvalue of the semisimple matrix B does 
not exceed 1 for a flat homogeneous Egorov coordinate system, then the signature of the given 
metric g a/3 (the difference between the numbers of positive and negative diagonal coefficients 
in a diagonal form of the metric) is minimal. (It is equal to for N even and to ±1 for N 
odd.) 

4. The third Hamiltonian structure of the Egorov homogeneous systems 

As was noted in Sec. 2, our physical examples of hydrodynamic type systems have two 
local Hamiltonian structures owing to the homogeneity and the invariance with respect to 
Galilean transformations. 

It is known from linear algebra that a skew-symmetric matrix B has pairs of eigenvalues 
differing in sign. In what follows, we sometimes do not distinguish between eigenvalues c a 
of B differing in sign. 

Let us expand the square of B as follows: 



ik 



5>' 

m^i,k 
m^i,k 



r i ){r m_ r k )pmip 



mk 



m\2, 



PmiPmk + W* + (r k ) 2 d k (3 kt + r^ lk + r k (3 kl 



- (r l + r K 



m^i,k 



+ r A r k 



r k )(P kl -P tk )- 



(39) 



^ PmiPmk + diPik + dkfiki 
■rriy^iyk 

Theorem 8. For even N = 2n, a homogeneous Egorov diagonal hydrodynamic type 
system (2) with flat metric g 1 ^ and second flat metric g%^ = r l g 1 ^ has a third local Hamil- 
tonian structure with metric g 1 ^ = (r % ) 2 g % ^ if and only if the matrix B is nonsingular and 
semisimple and has only one pair of eigenvalues ±c. 

Proof. In this case, B 2 = c 2 E, where E is the identity matrix, since an application of 
B 2 to any of its eigenvectors ip gives c 2 ip \ a \ On the other hand, it follows from (39) that in 
the presence of two local Hamiltonian structures determined by metrics g l fo and g 1 ^ = r l g 1 ^ 
(cf. (10) with /ij = r l ), one has 

[B\ k = £ (r m ) 2 P mi P mk + <?)%P* + {r k ) 2 8 k (3 kl + r% k + r k (3 kt = 

m^i,k 

provided that i ^ k, which precisely implies the existence of a third local Hamiltonian 
structure (cf. (10) with ji i ■ 



J\2\ 



Conversely, if condition (10) with n i 



J\2 



holds, then the matrix B 2 is diagonal. 



Since system (2) does not split, it follows that f3 ik ^ for % ^ k. By (17), it is obvious 
that none of the components of any homogeneous annihilator ip^ vanishes identically. Since 
B 2 vf)^ = c^ipf*^ , we see that the diagonal matrix B 2 has identical diagonal entries c 2 a . One 
can readily prove that the matrix B is semisimple by applying B 2 to root vectors; namely, 

if Bv = c a v + ip , then B v = c a v + 2c a ip ^ c a v. | 

Theorem 9. For odd N — 2n + 1, a homogeneous Egorov diagonal hydrodynamic type 
system (2) with flat metric g l fa and second flat metric g 1 ^ = r l g l fo has a third nonlocal 
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Hamiltonian structure with constant curvature metric g 1 ^ = (r l ) 2 g 1 ^ if and only if the 
matrix B is semisimple and has a simple eigenvalue C( ) = and a single pair of nonzero 
eigenvalues C(±&) = ±c (k = 1, . . . ,n) and the metric gJL itself is homogeneous of degree 
and depends only on the differences of the Riemann invariants r % . 

Proof. In this case, [-B 2 ]^- = c 2 5ij — (c 2 / g 00 )^^ ip^ , where ipf^ is a basis element in 
the kernel of the matrix B, Ripf^ = Bipf^ = 0, and g 00 = £\ ipf^ip^ is the corresponding 
element of the metric in flat coordinates. One can readily verify this identity by applying B 2 
to ipf*^ with regard to the form (38) of the metric in flat coordinates. Using the expansion (39) 
once more, we obtain 

E ( rm ) 2 PmiP mk + (rrdiP ik + (r k ) 2 d k (3 kt + r% k + r k /3 ki = 4fH®H<?\ 

since ipf^ = 2H^ (cf. (27)). This just means that there exists a third (this time, nonlocal) 
Hamiltonian structure determined by the metric gfy = (r l ) 2 g l fa of constant curvature c, 
since the expressions (13) computed for g l fo read 

-^-\B\ k = ^- h H^Hf\ i^k. 

r % r k L ltK r i r k % « ' ' 

The existence of the corresponding Hamiltonian for this nonlocal Hamiltonian structure 
readily follows from the results in [5]. 

Conversely, if an Egorov diagonal system with two flat metrics g l As and g l {L = r l g 1 ^ 
and with homogeneous functions /3 ik has a third nonlocal Hamiltonian structure with con- 
stant curvature metric g%s = (r*) 2 ^, then [B 2 ]^/ (r' l r k ) = di(r)Sij — cH i Hj/(r' l r k ), where 

Hi(r) = 1/ V g%} and the diagonal entries d^r) are some functions. Since the coefficients 
f3 ik are homogeneous and invariant with respect to the shift r l — > r l + const , it follows that 
all products HiHj are homogeneous of degree and invariant with respect to the shift. 
Consequently, expressing the functions Hi themselves via these products, we find that the 
metric gj[s is homogeneous of degree and depends only on the differences of the Riemann 

invariants r l . Applying the matrix B 2 to any eigenvector or root vector, we readily find that 
all nonzero eigenvalues are the same and that there are no root vectors. | 

Theorem 10. A homogeneous Egorov diagonal hydrodynamic type system (2) with flat 
metric g l fo and second flat metric g l L> = r l g 1 ^ has a third nonlocal Hamiltonian operator of 
the general form (11) with metric g l fo = (r l ) 2 g l fa provided that the matrix B is semisimple. 

Proof. In this general case, one has 

s 

where the annihilators ipf^ are basis elements of the kernel of the matrix B, the ipf*^ are 
its eigenvectors, and the ip[~ a ^ are its eigenvectors with opposite eigenvalues arranged in 
a way such that the metric (38) is diagonal for ipf^ and antidiagonal for ijj\ a \ijj\ a \ i.e., 
Y^i V^Vi"^ — 1; nere we nave chosen one of the eigenvalues c\ as the "main" diagonal entry 
of the matrix B 2 . This relation, just as in the preceding cases, can readily be verified by a 
straightforward substitution of the eigenvectors of B . Thus, we can use (27) and write out 
the corresponding expansion (13) guaranteeing the existence of the Hamiltonian operator 

ll 



(11). If B has the eigenvalues ±1/2 (see Sees. 5 and 6 below), one should take c\ = —1/2 
to express the desired annihilators ipf*^ via the corresponding h\°^ with the help of (27). | 
Note that Hf ] = } /2 can be found explicitly without quadratures even in the case of 
an arbitrary metric with homogeneous (3 ik . For example, for N = 3 we obtain ipf^ = (3 k j — 
r k (3 jk , and for A^ = 5 one has ip\ 0) = a jk a im + a jX a mk + a jm a kh where a jk = r 3 (3 kj - r k (3 jk . 
If the metric has the Egorov property, then we obtain ipf^ = (r J — r k )f3j k for iV = 3 and 
^(o) = {rj _ r k ){r i _ r m )(5 . k(3im + {r j _ r i)( r m _ r k)p jl( ] mk + (r i - r m )(r k - r l )(5 jm (5 kl for 

N = 5. In a similar way, one writes out ipf^ for TV = In — 1 > 5 . 

In the subsequent sections, we show how to establish the existence of second and third 
Hamiltonian structures for hydrodynamic type systems that are not written in Riemann 
invariants. 

5. Averaged iV-Phase Solutions of the Korteweg— de Vries Equation 
and the Nonlinear Schrodinger Equation 

Af -phase solutions of the KdV equation were averaged by the Whitham method in [22]. 
However, the Hamiltonian formalism for the averaged equations was developed later in 
[2,15,1]. The hydrodynamic type systems obtained by averaging inherit Galilean invari- 
ance, Hamiltonian property, and homogeneity. Moreover, since the KdV equation is Galilean 
invariant and hence has a pair of local Hamiltonian structures and a pair of conservation 
laws (see Theorem 1, Eq. (21)) 

d t u = d x [u 2 + e 2 u xx \, d t [u 2 + e 2 u xx ] = d x [^u 3 - 3s 2 u 2 x + 2e 2 (u 2 ) xx + e\ xxxx \, 

it follows that after the averaging on an iV-phase solution the resulting (2iV+ l)-component 
hydrodynamic type system is also Galilean invariant and has an Egorov metric (by The- 
orem 1) that is homogeneous and depends only on the differences of Riemann invariants. 
Hence, it has a pair of local Hamiltonian structures. The annihilators of the first local 
Hamiltonian structure of the averaged KdV equation comprise the averaged annihilator of 
the first local Hamiltonian structure of the KdV equation itself as well as iV wave numbers 
(corresponding to the iV phases 9i = kiX — oo{t of the quasiperiodic solution) and iV corre- 
sponding partial derivatives of the averaged Lagrangian with respect to these wave numbers. 
A straightforward computation shows that the homogeneity degree of the functions ipf^ 
corresponding to the averaged annihilator of the first local Hamiltonian structure is 0, the 
homogeneity degree of the functions ip[ k \ k — 1, . . . , N, corresponding to the wave numbers 
is —1/2, and the homogeneity degree of the functions ip'l N+k \ k = 1, . . . , N, corresponding 
to the partial derivatives of the averaged Lagrangian with respect to these wave numbers 
is 1/2. Thus (by Theorem 9) the averaged KdV equations have a third nonlocal Hamiltonian 
structure associated with a constant curvature metric (see [25]). (The existence of the first 
two local Hamiltonian structures was established earlier in [2].) 

Multiphase solutions of the nonlinear Schrodinger equation were averaged in the diploma 
paper of the first author. Later, this result was published in [6] . It was obtained independently 
in [21]. It was proved that averaged A^-phase solutions of the NLS equation have three 
local Hamiltonian structures and a fourth, nonlocal Hamiltonian structure defined by a 
differential-geometric Poisson bracket with constant curvature metric. The canonical form 
of the first three Hamiltonian structures in flat coordinates was given in [1] and [15]. The 
nonlinear Schrodinger equation written out in the Hasimoto form is Galilean invariant and 
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has three local Hamiltonian structures and the pair of conservation laws 



d t \u\ 2 = id x [uuz 



uu. 



.], id t [uu x - uu x ] = d x [\u\ 4 + A\u x \ 2 - [\u\ 2 ) xx ]. 



Reproducing the above argument for the KdV equation in the case of the nonlinear Schr6dingcr| 
equation, one can also show that these properties are inherited under iV-phase averaging 
(see [8]), which is a good illustration of Theorems 1 and 8. 

6. The Benney— Zakharov and Yajima— Oikawa— Mel'nikov Systems 

(the Dispersionless Limit) 



The Benney system (Zakharov's reduction; see [16,4]) 



w* = d x 



iu 



k\2 



N 



+ E 



171=1 



r/l = d x {u k r] k ) 



(40) 



is Galilean invariant (u k — > u k + c, x — > x — ct) and has a pair of conservation laws of the 
form (see Theorem 1, Eq. (21)) 

N 



( N 



N 



m=l 



^ m=l ' 

N , N 

Emv + 2 Ef 

m=l ^ m=l 



and the first local Hamiltonian structure 



1L 



a, 



6H_ 



Vt 



a 



determined by the following homogeneous Egorov metric depending on the differences of 
Riemann invariants (in the Riemann invariants \ k , k — 1, . . . , 2N; see [13]): 



0* 



E 



^ (//* + M m ) 3 ' 



where the characteristic velocities // and the Riemann invariants A* can be found from the 
system (see [23]) 

1] m ^ r^m 



E 



^ u % + u m ' ^ (u i + u m ) 2 

=1 ' m=l vl ' 

Comparing the homogeneity degrees of the flat coordinates (u k ,i] k ) and the metric g n , we 

(k) 

see that the homogeneity degree of the functions ip\ , k — 1, . . . , N , corresponding to the 
first half of flat coordinates u k is —1/2 and the homogeneity degree of the functions ip\ N+k \ 
k = 1, . . . , N, corresponding to the second half of flat coordinates i] k is 1/2. It follows by 
Theorems 1 and 8 that the Benney-Zakharov solution (40) has also second and third local 
Hamiltonian structures (see [10,9]). 

In the dispersionless limit, the Yajima-Oikawa system (known as the long-short reso- 
nance) generalized to the iV-component case (Mel'nikov's system I; see [26]) has the form 



d x 



k\2 



(u k ) 



+ w 



r] k = d x (u k r] k ), w t = d x 

13 



N 



E"' 

m=l 



It is also Galilean invariant (w — * w + c, x — > x — ct) and has a pair of conservation laws of 
the form 



/ N \ / N \ r N 

d t w = dJ^rrY d t [Y^r 1 m \=d x 



U m T] m 



■ m=l ' N m=l 

prescribed in Theorem 1 and the first local Hamiltonian structure 

k . 5H k a 5H 5H 

determined by the following homogeneous Egorov metric depending on the differences of 
Riemann invariants (in the Riemann invariants X h , k = 1, , . . . , 2N +1): 

Si 1 o \ / 



g n = i - 2 



^ (/ii + M™) 3 ' 
m=l v/ ( ' 

where the characteristic velocities fi { and the Riemann invariants A« can be found from the 

system 

2 N m N m 

m=l m=l v ' 

Comparing the homogeneity degrees of the flat coordinates (w,u h ,r] k ) and the metric g n we 
see that the homogeneity degree of the function ijjf^ corresponding to the flat coordinate w 

(k) 

is zero, the homogeneity degree of the functions ip\ , k — 1, . . . , N, corresponding to the 
first half of flat coordinates u k is —1/2, and the homogeneity degree of the functions ip\ N+k \ 
k = 1, . . . , N, corresponding to the second half of flat coordinates r] k is 1/2. It follows by 
Theorems 1 and 9 that the Yajima-Oikawa-Mel'nikov system has also second (local) and 
third (nonlocal, with constant curvature metric) Hamiltonian structures. 
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